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We give a general thermodynamic analyzis of the behaviour of the chemical potential of electrons in 
metals at a second order phase transition, including in our analysis the effect of long range Coulomb 
$Ih ' forces. It is shown, that this chemical potential can have a kink at Tc, both for fixed sample volume 

^ ' and fixed external pressure. The Coulomb term transfers the changes in chemical potential of the 

electrons into an experimentally observable shift of the surface potential if the sample is electrically 
CO ' connected to a ground potential. 
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i I. INTRODUCTION 

. The behaviour of the chemical potential (/i*^) of electrons at the superconducting phase transition was discussed 
recently in Ref. , and corresponding measurements were carried out in Ref. Q . It was shown theoretically, that 
in the BCS model /i^ is given by the expression /x*^ = — or, if the density of states near Ep is energy dependent 
^ /i"^ = — '^P/''-^ A(r)^), so that it has a kink at Tc- General consideration gives for the kink the expression 

o ■ 
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U ■ A{d^i''/dT)y _ dlnT, 



(1) 



' where is the number of electrons and ACy is the jump in specific heat. This conclusion was confirmed experimen- 
5h ] tally [Q, where it was shown that n'^ in YBa2Cu307_a; does indeed have a kink at Tc. This result is of an essentially 
. . . ■ thermodynamic nature and is actually independent of the specific nature of the phase transition, i.e. it is valid for 
other second order phase transitions as well, not only for superconductivity. 

However, in Ref. |^ relation Eq.^ was derived for a system at fixed volume, whereas the standard experiments are 
carried out at constant ambient pressure. For a system at constant presssure we should use the Landau expansion not 
for the Helmholz free energy F{y/r,N,ip), but the Gibbs free energy G{p, NjT^ip), where is an order parameter. 
Hence the well-known treatment [|| of phase transitions seems at first sight to rule out the possibility of obtaining a 
kink in fi: As the chemical potential is just this Gibbs free energy per particle, G — N/i, a kink in /i is the hall-mark 
of a first order phase transition hence it seems paradoxal that a kink in the chemical potential was observed 
experimentally in a second order superconducting phase transition. 

The situation, however, is not so simple, and the solution goes down to the definition of the chemical potential. The 
point is that, when dealing with metals, one usually discusses not the total chemical potential /i determined as a 
Gibbs free energy per elementary cell or per mole of the substance, but the chemical potential of the electrons /i"^ 
(which at T = is equal to the Fermi energy). One may say that the total chemical potential /i is the Gibbs free 
energy for adding to the system an extra unit cell, with all its nuclei and electrons, whereas /i"^ is the change in G 
when we change only the concentration of electrons. Thus, generally speaking, the situation here is similar to the 
situation in a mixtures, solutions or a complex chemical compound [^. In that case we have to introduce chemical 
potentials for each component 
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and whereas G itself, and the corresponding total chemical potential fj, = dG/dN = G/N (where N is the number 
of formula units for fixed concentration of components) should have no kink at Tc (this is actually the definition of a 
second order phase transition), partial chemical potentials /i^ may well behave differently, as has been discussed for 
the high ceramics by Burns [Q . The situation with n'^ is still somewhat more complicated because, if we want to 
discuss it separately, i. e. if we want to look at the changes in the system when we change the electron concentration, 
we have also to take into account long range Coulomb forces, which usually guarantee electro-neutrality of the system. 
We will show below that, all these factors taken care of, the chemical potential of electrons /i*^ as measured e.g. in 
Ref M should indeed have a kink at Tc, even if in the experiment not the volume but the pressure is fixed externally. 



II. SINGLE COMPONENT SYSTEM. 



First of all we derive Eq. |^ using a slightly different method than that used in ||^ , and show that indeed even for 
a one-component system the chemical potential at constant volume has a kink at the second order phase transition, 
in contrast to the situation at fixed pressure. At fixed pressure the Gibbs free energy G{p,T,N) = Nfi{p,T), and 
/i(p, T) as well as the first derivatives of /i are continuous at the transition. The temperature derivative of n at fixed 
pressure follows by making the transformation of variables T) = fi(p(V,T),T), so that 
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We furthermore identify {dfi/dp)rp as V/N. Only the second term of this expression gives rise to a discontinuity at 
the phase transition, for which we can use one of the Ehrenfest relations Q 

A {dp/dT)y ^ d\nT, 

ACv dN ^ ' 

As for a single component system NdTc/dN — —VdTc/dV, we can now rewrite the Ehrenfest relation in the form 



A{dfi/dT)y _d\nTc 



ACv dN 
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which is identical to Eq.n^. This treatment resolves the apparent contradiction between microscopic treatments, which 
give a finite jump in dji/dT at Tc, and general arguments according to which ^ has to be smooth at a second order 
phase transition: As the microscopic treatments are carried out for a given electron concentration, and neglecting 
temperature dependencies of the volume, if the system is held at a constant external pressure. As we will show 
in the next section, the partial chemical potential of electrons may have a kink, even if such thermal changes of 
the volume are taken into account. Thus it is instructive to look for its behaviour first at fixed volume, such as 
is usually obtained from a microscopic description, and than make the Legendre transformation of the free energy 
G{p,N,T) = F(V,N,T) +pV, to discuss the effects of fixing external pressure instead of volume. Important here is 
the role of the long range Coulomb interaction 



III. FIXED VOLUME 



The thermodynamical state of the solid as follows from these conditions is described by the Helmholtz free energy 
F{V,T, N, Ne), where N is the number of elementary cells, and iVe is the number of electrons. Equivalently we can 
use the electron concentration x = N^/N, so that we can write the free energy in the form F{V,T,N,x). The Gibbs 
free energy per unit cell is fJ.{V) = (f^)^^.^. and the chemical potential of the electrons is ^J.'^{V) — N^^ {^)vt n' 
If the particles would be neutral, the free energy would depend on the volume a,s F — Vf{n,x,T) where n — N/V. 
However if we want to discuss phenomena where the number of electrons is allowed to change, we should also add the 
Coulomb (charging) term, which has a different dependence on the size of the system. We use for the corresponding 

term the simple expression SFc — ^ ^2V^^°'' (Coulomb energy of a charged sphere of radius R = V^^^ and charge 
eN{x — xo) with xq the equilibrium number of electrons per unit cell). We will see later, that the precise shape of the 
sample, and therefore the exact value of the prefactor in 6Fc is unimportant. Hence the Helmholtz free energy is 
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F{V,N,x,T) = Vf{n,x,T) + ^e^{x-xo)'^N^V-'^/^ (6) 
For systems where the external volume is fixed, the chemical potential of the electrons is 

H%n,x,T,cl^) = ^|||^^^^^ = i|£+e0 (7) 

where in the last line we introduced the charging potential (f> = e{x — xo)n^^^N^^^ ■ If the sample is externally 
grounded, the electrons are in equilibrium with the ground potential which fixes fi'^ at the value of the external bath. 
The electrons can now flow freely in and out of the sample, so that the number of electrons follows from solving the 
above expression for x. From this we see that the charge per unit cell x — xq is proportional to N~'^/'^. Note that the 
total chemical potential has a different form: 

We see that, due to the fact that (p is finite, the last (Coulomb) term vanishes in the thermodynamic limit. Clearly, 
as different partial derivatives are taken of the functional /(n, x, T), also a different behaviour at the phase transition 
occurs for /i and /i*^. Both of them, however, would have a kink at Tc: Using the continuity of entropy along the curve 
Tc{x) in the former, and Tc{n) in the latter case (as in Ref. [^), one obtains Eq. |l| for the kink in d^'^/dT, while 

Hence the kink in dfj,/dT is proportional to the derivative of Tc with respect to the density of unit cells, which may 
also include a possible redistribution of electrons between reservoirs upon changing the lattice constant. On the other 
hand A(dfi'^/dT) is proportional to the dT^/dx, where x is the number of electrons per unit cell. As was treated 
in Ref. [0|, if the solid contains two charge reservoirs, one of which is 'active' in the phase transition, a reduction of 
A (dfi'^ / dT) occurs if a charge re-distribution between the reservoirs takes place at the phase transition. This also 
follows from the present thermodynamic analysis: dTc/dx will be reduced, if part of the electrons moves into an 
'inactive' reservoir upon varying x. 



IV. FIXED PRESSURE 



First of all we again obtain the general thermodynamic expression for the behaviour of /x*^ and {dfi'^ / dT)^ initially 
ignoring charging effects. Following |^ we differentiate the expression A (dG/dT) = (continuity of entropy) along 
the curve Tc{x) and find 
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Thus we see, that indeed if we could treat the electron concentration x as an independent parameter, the chemical 
potential of electrons fi*^ would have a kink at the second order phase transition, even at fixed pressure. Expression 
pj] has the same structure as Eq. ||, with the natural change of variables. 

We can take into account charging effects, by going from F{V, N, x,T) to the Gibbs free energy G ^ F +pV, where 
p — —dF/dV, and calculate fj,'^{p, N,x,T) = dG/dNe and uip,N,x,T) = dG/dNe- As a result of the Legendre 
transformation, these have the same functional form as Eqs. and |[ However, these expressions are in this case a 
fmiction of p, as the dependence of the density n on the external pressure has to be solved from 



p{n,x,T) = -f{n,x,T) + + U^n^^N-y^ 



(12) 



As before 0, defined as e{x — xo)it}^^ N'^/^ , is the charging potential. Note that the Coulomb term introduces an 
additional dependence on the total number of unit cells N , due to its long range nature. Hence in principle we should 
consider n{p, T, x, N). However, this contribution vanishes for N ^ oo. 
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_^ ^2^-1/3^-2/3 



(13) 



As fi has to be evaluated at the minimum with respect to the order parameter, it is automatically ensured that 
it has no kink at the phase transition, provided that a Landau expansion can be made. We also notice, that the 
charging effects do not affect our expression for the kink in fj,^ (Eq. [Tl]) , as the charging term is of vanishing order 
for N oo: If the sample is electrically isolated N and x are fixed. As (p — e{x — xo)n^^^N^^^ , the only thermal 
variations enter through the volume changes (changes of n). Hence, if the sample is electrically charged, there are 
temperature dependent changes to fi"^ due to changes in lattice constant at fixed external pressure, that should be 
taken into account. If, on the other hand the sample is electrically grounded, ^'^ is fixed externally. In this case ecj) 
is exactly equal to /i"^ (but with opposite sign) of the electrically isolated uncharged sample, and can be measured 
exp er iment ally. 



V. EXAMPLE 

For a solid undergoing a second order phase transition, we may write f{n, x, T, ip) 
we make the following free energy expansion 



fo + U, where for /o and 



Fo(7V,t/,T,V) = (Mo + /iv)^ + 5^oK)(^-l)' + c(V-K))|V^P + i</)2yi/3 
/xV,(n,T,^) = a(T -r,)|i/;|2 + ||^|4 



(14) 



It has to be understood here, that the parameters a, b, and are independent of the volume V, and coupling between 
the strain field and the superconducting order parameter is introduced in an ad-hoc manner through the constant c. 
The expression for the pressure as occurs in Eq. n3 is 
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6^ 



(15) 



from which we solve the volume V^ = Vb{l — (p-t- c\ip\'^) / Bq^ . We are now ready to calculate the Gibbs free energy, 
and from it, with the definition of /i, the chemical potential 
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Here no = N/Vq represents the molar density of the solid if the order parameter is zero. So in the first place we 
notice, that in spite of the coupling to the lattice, the transition remains of second order. Only if the coupling constant 
c exceeds the critical value (bnoBoY^^ , it becomes of first order. In the second place we notice that, in the case where 
the transition is still of second order, Tc is shifted with an amount and direction which depends on c and the value of 
the externalUy applied pressure. The Gibbs free energy per particle calculated at the minimum with respect to i/'l^ is 
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Clearly the total chemical potential of the system has a smooth behaviour at the phase transition without a kink, 
unless the coupling to the lattice is sufficiently strong to make the transition first order. 

Let us now consider the electronic sub-system. As we are considering a solid, the number of unit cells, N, is externally 
fixed, and the corresponding chemical potential need not be in equilibrium with an external bath. The situation is 
completely different for the electrons however, as these can move in and out of the solid. Using again the definition 
X = Ne/N for the number of electrons per unit cell, we can calculate the corresponding chemical potential by 
differentiating the total Gibbs free energy {^iN, where /z is the chemical potential as calculated above) with respect 
to Ne- As we have seen in the preceeding section, the long range Coulomb forces can be included, and the correct 



expression becomes fJ-'^{p) = 



ecj), where we have to assume now that /i not only depends on p and T, but 



also on X. Indeed most of the properties of a solid depend strongly on the number of charge carriers per unit cell, for 
example by influencing the strength of the chemical bond between neighbouring unit cells, or by having an effect on 
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the superconducting transition temperature. For our discussion the latter dependence is the most important one, as 
we are interested in the behaviour near the superconducting transition. Let us indeed assume that /i^ is derived from 
a microscopic theory, which also predicts that dTc/dx ^ 0. If the electronic subsystem is brought into equilibrium 
with an external bath by connecting it with a current lead, the charge on the sample is such, that e(j) compensates for 
the difference. Hence /i'^(p) = /i'^^* and the voltage on the sample is 

e(j) = ii - ? 1— 18 

h dx 

This voltage has a kink at Tc, which can be determined experimentally by measuring the workfunction of the 
sample. If the lattice is sufficiently soft {Bq < c^/nob) the transition becomes of first order. The above expression 
diverges at the point where Bq — c^/nob and is no longer valid in this limit. If the material is hard {Bq 3> c^/nob) 
and the external pressure small {p ^ noaBo/c) Eq. HS^ is just the result which we already obtained at fixed sample 
volume. 



VI. CONCLUSIONS 



The seemingly paradoxal result of microscopic theories of superconductivity, that the chemical potential may have 
a kink at a second order phase transition, is resolved. There are two ingredients in the resolution of the paradox. 
The first is, that the microscopic treatment is always carried out at a given concentration of electrons, i.e. at fixed 
volume or density. In contrast to fixed pressure, there is absolutely no general rule forbidding the kink in /x at fixed 
volume; general thermodynamic results ||] confirm that. More interesting is the second part of the story. The chemical 
potential of the electrons is in general different from the total chemical potential and thus can have a kink, even at 
fixed pressure. This stems from the fact that electrons are free to move in and out of a solid, thus maintaining electrical 
equilibrium with the environment, whereas the ions in a solid are immobile. As a result the chemical potential of 
electrons and the Gibbs free energy per mole of the ions have different dependencies on temperature, pressure etc., 
and are only coupled through the long range Coulomb forces. Using scaling arguments the Coulomb charging energy 
is shown to be of vanishing order in the Gibbs free energy per unit cell, the density of the solid, Tc and the order 
parameter. At the same time the long range Coulomb forces on the one hand keep the charge carrier density fixed, 
while on the other hand they transform changes of the chemical potential of the electrons into equally large and 
measurable changes of the workfunction, if the solid is in electrical equilibrium with its environment. 
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